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Abstract 

Gene silencing and RNA interference are major cellular processes that control gene expression via the cleavage of target 
mRNA. Eukaryotic translation initiation factor 2C2 (EIF2C2, Argonaute protein 2, Ago2) is considered to be the major player 
of RNAi as it is the core component of RISC complexes. While a considerable amount of research has focused on RNA 
interference and its associated mechanisms, the nature and mechanisms of nucleotide recognition by the PAZ domain of 
EIF2C2/Ago2 have not yet been characterized. Here, we demonstrate that the EIF2C2/Ago2 PAZ domain has an inherent 
lack of binding to adenine nucleotides, a feature that highlights the poor binding of 3'-adenylated RNAs with the PAZ 
domain as well as the selective high trimming of the 3'-ends of miRNA containing adenine nucleotides. We further show 
that the PAZ domain selectively binds all ribonucleotides (except adenosine), whereas it poorly recognizes 
deoxyribonucleotides. In this context, the modification of dTMP to its ribonucleotide analogue gave a drastic improvement 
of binding enthalpy and, hence, binding affinity. Additionally, higher in vivo gene silencing efficacy was correlated with the 
stronger PAZ domain binders. These findings provide new insights into the nature of the interactions of the EiF2C2/Ago2 
PAZ domain. 
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Introduction 

Gene silencing and RNA interference (RNAi) are cellular 
processes that control gene expression via the cleavage of target 
mRNA by either a cellular or synthesized small RNA [1,2]. A 
small interfering RNA (siRNA) is loaded into an RNA-induced 
silencing complex (RISC) that guides the RISC complex to bind 
with the complementary mRNA sequence followed by endonu- 
cleolytic cleavage of the target mRNA, thus preventing gene 
expression [3-5] . RNAi technology has been found to be a useful 
strategy in the fight against cancers [6-8], viral infections [9-12], 
asthma [13,14], rheumatoid arthritis [15], and other diseases 
[16,17]. Eukaryotic translation initiation factor 2C2 (EIF2C2, 
Argonaute protein 2, Ago2) is an Argonaute (Ago) protein and 
considered to be the major player of RNAi. 

Human Ago proteins are classified into the Ago and PIWI 
subfamilies. The Ago subfamily consists of four closely related 
proteins, namely, Agol (EIF2C1), Ago2 (EIF2C2), Ago3 
(EIF2C3), and Ago4 (EIF2C4), while the PIWI subfamily includes 
HIWI, HILI, PIWIL3, and PIWIL4 [18]. Ago proteins have a 
conserved structure that is composed of two binding domains 
(PAZ and MID) and one catalytic (PIWI) domain [19]. The PIWI 



domain contains an RNase H fold and is considered as the 
catalytic RNA cleavage center of Ago proteins. The 5'-end of 
siRNA is anchored in the cavity of the MID domain, while the 
PAZ domain binds to the 3 '-double nucleotide overhang of siRNA 
[20-24]. Interestingly, the function of the PAZ domain extends 
beyond a simple binding event to include significant effects on 
RNA loading and the sheer (endonuclease) activity of the RISC 
complex. PAZ-disrupted mutants of Ago proteins are unable to 
unwind or eject the passenger strand and prevent the formation of 
active RISC complexes [25]. Furthermore, chemical modifications 
of nucleotides that bind with the PAZ domain affect the overall 
RNAi process [26-28]. Within the Ago proteins of humans, 
EIF2C2/Ago2 is the only member with recognized slicer activity 
and is considered as the catalytic engine of the RNAi process. 
Furthermore, cells lacking EIF2C2/Ago2 are unable to process 
siRNA [29]. 

The characteristics and binding parameters of die PAZ domain 
(named after the three proteins that were initially found to contain 
this domain, such as PIWI, Ago, and ZwUle) from Ago proteins 
have been characterized in Dro.sophila melanogaster Ago 1 and Ago2 
[22-24,30], human EIF2Cl/Agol [21], human HIWI [31], HILI, 
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PIWIL3, and PIWIL4 [32], and mouse PIWILl [33]. However, 
the exact single nucleotide specificities of the PAZ domain are 
not well understood. The interaction of the 3 '-end of RNA with 
the PAZ domain shows interesting molecular dynamics. During 
RNAi, the 3'-end of RNA toggles between binding and release 
from the PAZ domain in a cyclic manner. This process was 
demonstrated to be essential for proper RNAi and cleavage of 
the target RNA. In this context, we recently evaluated the 
impact of the strength of binding with the PAZ domain on the 
overall RNAi process using a computational approach [34,35]. A 
small interaction surface and lower free binding energy were 
among the markers for favorable RNAi. To this end, there are 
still unknown factors in the recognition of RNA by the PAZ 
domain and several aspects of the modified 3 '-end of siRNA 
have yet to be determined. Of special interest, the nature of 
nucleotides at the 3'-dangling end of siRNA is not well 
understood. In this study, we asked several questions to evaluate 
the exact molecular aspects of the PAZ domain. First, which of 
the known nucleosides can bind with the PAZ domain? Second, 
what are the determinants of nucleotide recognition by the PAZ 
domain? Lastly, what is the relationship between strong and 
weak binders and their in vivo RNAi efficacy? Here, we 
demonstrate that the EIF2C2/Ago2 PAZ domain has a selective 
binding affinity for ribonucleotides over deoxyribonucleotides. 
Furthermore, there was an inherent bias against adenine 
nucleotides. Within all ribonucleotides (rA, rG, rC, and rU), 
unexpectedly, rA showed very low affinity compared with the 
other ribonucleotides, which constitute the highest affinity 
ligands. We further show that the EIF2C2/Ago2 PAZ domain 
functions as an induced fit cavity and adapts to different 
nucleotides with no or minimal changes in conformation. 
Finally, we demonstrate that the higher binding nucleotides 
were generally associated with higher in vivo RNAi activity. 



Materials and Methods 

Design and synthesis of siRNAs 

The used siRNAs were designed to suppress expression oiRmilla 
luciferase. The 3' overhang contained a single base extension of 
natural ribonucleotides or deoxynucleotides (Table 1). The siRNAs 
1-8 (oligos no. 1-32) were purchased from Hokkaido system 
science co. (Hokkaido, Japan). Oligos no. 17 and 18 were 
synthesized witli DNA/RNA synthesizer by using phosphoramidite 
method. Deprotection of bases and phosphates was carried out 
night in NH40H/EtOH (3:1) solution. 2'-TBDMS groups were 
removed by 1 M tetrabutylammonium fluoride in THF for 
1 2 hottrs. The reaction was stopped by 0. 1 M TEAA buflTer pH 7 
and desalted on Sep-Pack C18 cartiridge. Deprotected oligos were 
purified by 20% PAGE containing 7 M urea. The concentration of 
oligos was measured as OD units at 260 nm. The extension 
coefficient of the oligos was calculated from the mononucleotides 
and dinucleotides according to the nearest-neighbor approximation 
method. MALDI-TOF/MS analysis of oligos was performed to 
detect the mass of oligos with a time-of-flight mass spectrometer. 
siRNA9 was used in routine gel shift assays during purification of 
PAZ domain. 

Cell cultures and transfection 

HeLa cells were grown at 37°C under optimal humidity and 5% 
CO2 in Dulbecco's Modified Eagle's Medium (D-MEM) supple- 
mented with 10% bovine serum (BS). Twenty four hours before 
transfection, cells were diluted to 40000 cells/ml and aliquots of 
cells culture were transferred into 96-well plate (100 |j,l/well). Cells 
in each well were transfected with 35 |j.l of 20 ng psiCHECK-2 
vector, 1 nM or 10 nM of each siRNA (siRNAs 1-1 6), and 1.2 |ig 
of transfection reagent in Opti-MEM Reduced Serum Medium. 
Cells were incubated for 24 hours and subjected to analysis of 
gene silencing efficacy. Control wells were adopted by transfection 
of cells in a mediun lacking siRNAs. After 1 hour, 100 |ll MEM 



Table 1. Sequence of oligonucleotides used in tliis study. 



No. of SiRNA 


Oligo. No. 


Ollgo 


Sequence 


No. of ON 


siRNAl 


ON 1 


Renilla-rG AS 


5'-GUAGGAGUAGUGAAAGGCCG-3' 


20 




0N2 


Renilla-rG S 


5'-GGCCUUUCACUACUCCUACG-3' 


20 


siRNA2 


0N3 


Renilla-rU AS 


5'-GUAGGAGUAGUGAAAGGCCU-3' 


20 




0N4 


Renilla-rU S 


5'-GGCCUUUCACUACUCCUACU-3' 


20 


siRNA3 


0N5 


Renilla-rC AS 


5'-GUAGGAGUAGUGAAAGGCCC-3' 


20 




0N6 


Renilla-rC S 


5'-GGCCUUUCACUACUCCUACC-3' 


20 


siRNA4 


0N7 


Renilla-rA AS 


5'-GUAGGAGUAGUGAAAGGCCA-3' 


20 




0N8 


Renilla-rA S 


5'-GGCCUUUCACUACUCCUACA-3 ' 


20 


siRNAS 


0N9 


Renilla-dG AS 


5'-GUAGGAGUAGUGAAAGGCCdG-3' 


20 




ONIO 


Renilla-dG S 


5'-GGCCUUUCACUACUCCUACdG-3' 


20 


siRNA6 


ONll 


Renilla-dT AS 


5'-GUAGGAGUAGUGAAAGGCCT-3' 


20 




ONI 2 


Renilla-dT S 


5'-GGCCUUUCACUACUCCUACT-3' 


20 


siRNA7 


ONI 3 


Renilla-dC AS 


5'-GUAGGAGUAGUGAAAGGCCdC-3' 


20 




0N14 


Renilla-dC S 


5'-GGCCUUUCACUACUCCUACdC-3' 


20 


siRNAS 


0N15 


Renilla-dA AS 


5'-GUAGGAGUAGUGAAAGGCCdA-3' 


20 




0N16 


Renilla-dA S 


5'-GGCCUUUCACUACUCCUACdA-3' 


20 


siRNA9 


ONI 7 


Renilla-F-dT AS 


5'-F-GUAGGAGUAGUGAAAGGCCT-3' 


20 




0N18 


Renilla-F-dT S 


5'F-GGCCUUUCACUACUCCUACT-3' 


20 



doi:l 0.1 371 /joumal.pone.0094538.t001 
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contaimng 10% BS was added to each well and incubated for 
24 hours. 

Measurement of gene silencing 

The ability of siRNAs to prevent gene expression was 
determined by dual-luciferase assay using psiCHECK-2 vector, 
which contains Rmilla and firefly luciferase gene. The used siRNAs 
are targeting Renilla luciferase gene. Cell extracts were prepared in 
passive lysis buffer. The activities of Renilla and firefly luciferases in 
cell lysates were determined by a commercial Dual-Glo luciferase 
reporter assay kit. The signal of Renilla luciferase was compared 
with that of firefly luciferase from three different transfection 
experiments each contain three cultures for every oligo. 

Cloning, expression and purification of Ago2 PAZ domain 

The Ago2PAZ domain was cloned into the BamHI and Hindlll 

sites of pQE30 vector by using the standard cloning protocol. To 
generate PAZ domain, a fragment of Ago2 (residues 226-379) was 
amplified by using a pair of primers — forward (5'- ATAG- 
GATCCGCACAGCC AGTAATCGAGTTTGTTTG-3 ') and 
reverse (5'- GCGAAGCTTAATCTCTTCTTGCCGATCGG- 
GC-3'). The recombinant E. coli cells were grown in LB medium 
containing 50 Hg/ml ampiciHin at 37°C to the mid-log phase 
(-Deoo = 0.6). Induction of expression was carried out by the 
addition of IPTG (isopropyl P-D- 1 -thiogalactopyranoside) to a 
final concentration of 0.5 mM and cell growth was continued at 
20°C for 8 h. Cells were lysed in an extraction buffer (25 mM 
potassium phosphate buffer pH 7). After extraction, ethyleneimine 
polymer was added to a final concentration of 0.05% to precipitate 
nucleic acids. The hexahistidine-tagged PAZ domain was purified 
from the soluble cell extract by using talon metal affinity resin 
followed by gel filtration chromatography by using sephacryl 
S200HR. A final polishing step was performed by using Amicon 
Ultrafree-MC centrifugal filter unit (5 kDa membrane cut off 
value). 

Isothermal titration microcalorimetry 

Calorimetric experiments were carried out using VP-ITC. 
Before ITC experiments, PAZ domain was thoroughly dialyzed 
against ITC buffer composed of 25 mM potassium phosphate 
buffer pH 7. After dialysis, PAZ solution was centrifuged at 
15000 RPM for 20 min and the precipitate was discarded. 
Concentrated PAZ solution was kept at 4°C and diluted with 
the final dialysis solution immediately before use in ITC titrations 
to a final concentration of 30 |iM. Substrate solution was made 
from the final dialysis buffer of proteins to minimize artifacts due 
to subtle differences in buffer composition. The protein solution 
was thoroughly degased and loaded into the sample cell. The 
substrate solution was loaded into the syringe and used to titrate 
the protein solution by twenty sequential injections. Control 
experiments were performed by injecting the nucleotide into the 
dialysis buffer to determine the heats of dilution. The apparent 
heat change after each injection was integrated and corrected for 
the heat of dilution of the nucleotide. The data was fitted to a 
single-binding site model by non-linear regression analysis to yield 
the thermodynamic paramcti^rs K],, association constant; A//, 
enthalpy of binding; and n, the stoichiometry of binding. The 
affinity of the nucleotides to protein is given as the dissociation 
constant (Kj= \/K^. The binding entropy contributions were 
calculated from the equation AG= AH-TAS and AG = -RT In K, 
where Kis the association constant, Tis the absolute temperature, 
AS is the entropy change, and AG is Gibb's free energy change. 



Heat capacity measurement 

To address the association between structural and thermody- 
namic data, ITC titration experiments were conducted at four 
different temperatures between 8 and 3()°C and the thermody- 
namic parameters were determined. The change in heat capacity 
was determined from the temperature dependence of the enthalpy 
according to the following equation: ACp = 8AH/8AT(l) 

Mutual competitive displacement experiments 

Free PAZ domain was titrated wdth AMP, dAMP, dCMP, 
uridine or thymidine until reaching saturation. After finishing the 
first experiment, the ITC syringe was removed, thoroughly 
washed and loaded with UMP and a second titration experiment 
was carried out. All competitive binding experiments were carried 
out at 12°C. For simulation of competitive binding data, analysis 
model was used as previously described. 

Structure-thermodynamic calculations 

To test whether the binding of nucleotides with PAZ domain is 
associated with changes in the proteins accessible surface area, we 
analyzed several crystal structures by using VADAR web server 
(Volume, Area, Dihedral Angle Reporter). The ACp was calculated 
using the following equation: AC,, = 0.45 AAASAapolar-0.26 
AAASApolar (2), where AAASAapolar and AAASApolar are the 
changes in ASA of polar and apolar residues, respectively. 

From the output of data, we .selected and compared the total 
and fractional values of surface polar and nonpolar surface areas. 
The obtained data were used in empirical equations to analyze the 
change in surface area and its components after binding with the 
substrates. 

Synthesis of 5-methyl uridine monophosphate (rlMP) 

During the ITC assay, the binding afiinity of deoxythymidine 
nucleoside was surprisingly low. Furthermore, dTdT is a standard 
3'-terminal additive in many RNAi assays. Binding affinity studies 
revealed that the PAZ domain bears an inherently high affinity for 
ribonucleosides (except rA). Therefore, we designed a new ITC 
experiment to check the hypothesis that the PAZ domain 
preferentially binds ribonucleosides (except rA), regardless of the 
attached base. Moreover, this modification will give us an insight 
into the functional role of the 2' -OH of nucleotides. In this 
context, 5-methyl uridine was considered as a thymidine base 
attached to a ribosugar. Thus, a new comparison was carried out 
between deoxythymidine and ribothymidine nucleosides. 

5-Methyluridine-5'-0-monophosphate was prepared as follows. 
Distilled HgO (28 |tL, 1.56 |tmol), MeCN (2.38 mL), pyridine 
(215 ^iL), and POC13 (206 |J,1, 2.21 mmol) were added to a 20-mL 
round-bottom flask on ice, then 5-methyluridine (1) (100 mg, 
387 nmol) was added. The mixture was stirred at 0°C under argon 
gas for 5 h. The reaction was moved to room temperature and 
stirred for an additional 43 h. The reaction was checked by TLC 
(BuOH: H20: acetic acid =2:2:1), and quenched by iced H20 
(1 mL). The mixtur(; was concentrated and the residue was 
purified by ion exchange column chromatography (0-5% formic 
acid in H20), then the eluate was neutralized by NH40H (1 mL). 
The desalted supernatant solution was concentrated to give 5- 
methyluridine-5'-0-monophosphate (2) (98.4 mg, 291 jmiol, 
75%) as a brown solid. The resulting compound was checked by 
NMR (Figure SI). 

Statistical analysis 

The statistical analysis of RNAi potency was carried out using 
STATA version 12.1. A single-tailed one-way analysis of variance 
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Figure 1. Representative ITC profiles of the binding of nucleotides to the free PAZ domain. The top panels show the raw heat change 
data. The lower panel shows the integrated binding isotherms as a function of the molar ratio of the ligand to the enzyme. 
doi:1 0.1 371 /journal.pone.0094538.g001 



with Boiiferroni's multiple comparison test was conducted. The 
results of statistical analysis are provided in Tables SI and S2. 

Other methods 

Protein concentration was determined using a protein assay 
reagent with bovine serum albumin as the standard. The 
molecular mass of the recombinant enzyme was determined by 
SDS-PAGE (15% polyacrylamide). After electrophoresis, the gel 
was stained with Coomassie Brilliant Blue R-250. For Western 
blotting, purified protein was applied to 15% SDS-PAGE, 
transferred to a PVDF membrane, and detected by anti-His tag 
(10,000-fold dilution). The immunoprobed proteins were detected 
using an HRP conjugated substrate detection kit. 

Results 

Isothermal titration calorimetry 

ITC has largely became the method of choice in determining 
the binding affinity of two interacting partners [36-39] . Further- 
more, in one experiment we can get the full thermodynamic 
profile associated with binding. In this work, ITC was carried out 
for determining the binding parameters as weU as the mechanisms 
of interaction of different nucleosides with Ago2PAZ domain. The 
typical ITC titrations included addition of the protein in the cell 
and titrated with the ligand in a time controlled process. So that, if 
a binding event is likely to occur, the initial titrations wiU give high 
exothermic or endothermic peaks. As the titrations proceed, the 
binding sites will be gradually occupied by the ligand giving rise to 
smaller peaks. In the last injections, as the binding sites are 
completely occupied, the heat arising from injections are similar to 
the heat of background. Figure 1 shows the typical titration of 
Ago2PAZ domain with the nucleotides. The top panel shows the 
raw calorimetric data referring to the amount of heat produced 
following each injection. The bottom panel shows the integrated 



amount of heat generated per injection as a function of the molar 
ratio of ligand to protein. The obtained thermodynamic constants 
are given in Table 2. 

The isothermal titration calorimetry (ITC) assays included the 
ribonucleoside monophosphates UMP, GMP, CMP, and AMP, 
the deoxynucleosides monophosphates dTMP, dGMP, dCMP, 
and dAMP, and nucleosides such as uridine and thymidine. No 
binding activity was measured with thymidine or uridine (Fig. 1, 
Table 2), indicating the importance of nucleoside derivatives; for 
example, phosphate, in the recognition of substrates by the PAZ 
domain. Within ribonucleotides, binding activity was detectable 
with UMP, GMP, and CMP. Surprisingly, AMP binding was not 
measurable, indicating no or very weak binding of AMP with the 
PAZ domain. Within the deoxyribonucleotides, dGMP showed 
the most marked binding affinity. Furthermore, the nucleoside 
thymidine, which is used as a standard nucleotide at the 3'- 
overhang of siRNA in many RNAi experiments, showed very low 
binding activity. 

Thermodynamic signature. The thermodynamic signa- 
tures of UMP, GMP, CMP, TMP, and dGMP binding with the 
EIF2C2/Ago2 PAZ domain indicated different binding patterns 
among the ribonucleosides and deoxyribonucleosides (Figure S2). 
The thermodynamic signature of nucleosides included a negative 
AG, positive TAS, and negative AH for UMP, GMP, and CMP, 
and a positive AH for dTMP and dGMP. A negative AH value 
indicates favorable enthalpic conditions, while a positive value 
indicates unfavorable enthalpic conditions. Furthermore, a posi- 
tive TAS value indicates favorable entropic conditions. Favorable 
enthalpic and entropic terms are the most favorable thermody- 
namic profiles for tight binding. Therefore, the weak binding of 
dTMP and dGMP is favored by adverse enthalpic conditions due 
to their positive AH values. 

Substrate binding affinity. UMP showed the highest 
binding affinity (with Kj value of 15 |J.M) followed by GMP, 
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CMP and dGMP, respectively (Figure S3, Table 2). TMP showed 
at least 1 5-fold increase in the dissociation constant compared with 
UMP. 

Binding enthalpy and entropy. The interaction of ribonu- 
cleotides (except rA) was characterized by low negative enthalpy 
with A// value <1 kcal/mol (Table 2). In contrast, the deoxyr- 
ibonucleosides dTMP and dGMP showed a low positive enthalpy 
change. The positive AH value indicated unfavourable enthalpic 
conditions. 

Mutual competitive binding experiments 

In a mutual competitive binding experiment, we assumed that 
there was residual binding of the ribonucleotide AMP and the 
deoxyribonucleotides dAMP and dCMP or uridine nucleoside 
with the PAZ domain. Thus, two sequential experiments were 
performed. The first experiment included the titration of the first 
(low or no affinity) ligand followed by titration with UMP. If 
binding of AMP, dAJVIP, dCMP, or uridine is likely to occur, the 
binding of UMP could be minimized by the degree of binding of 
the other ligand. As a control, the same experiments were carried 
out in the presence of stronger binding ligands such as GMP. 
Figure 2 shows the results of the competitive binding experiments. 
Titration of the EIF2C2/Apo2 PAZ domain with UMP gave the 
characteristic figure shown as (UMP). Titration with UMP into a 
preformed complex of the PAZ domain with GMP (GMP-UMP) 
or CMP (CMP-UMP) showed very weak signal without measur- 
able thermal parameters, corresponding to a lack of UMP binding. 
Furthermore, in the presence of the medium affinity binding 
ligand dGMP (dGMP-UMP) or TMP (TMP-UMP), UMP showed 
a corresponding decreased access to the binding site as evidenced 
by lower exothermic peaks compared with the control (UMP). 
Finally, there was a slight decrease in the affinity of UMP in the 
presence of other ligands (dCMP, dAMP, and uridine). The 
calculated values of the weak binding ligands are given in Table 3. 
The estimated values of such low affinity ligands were 72- 
159 M while uridine was not expected to bind with the PAZ 
domain. 

UMP/GMP-induced heat capacity change (ACp) 

The measurement of heat capacity is an important parameter in 
determining the structural changes of the PAZ domain in response 
to its ligands. ACp can be determined by running the typical ITG 
experiment at different temperatures. For this purpose, the PAZ 
domain was titrated with UMP or GMP and the change in 
enthalpy was recorded at the different temperatures (Table 4). 

At lower temperatures, the binding of UMP or GMP was under 
enthalpic forces due to their negative AH values. Furthermore, 
there was a gradual decrease in enthalpic forces as the temperature 
increased. The enthalpy changes of UMP AH™''/AH^^'' and 
AH^^'' /AH^^wae 13.6 and 3, respectively. Thus, the enthalpy 
change decreased by approximately 1 3-fold when the temperature 
was raised from 8 to 30°C. The decrease in enthalpic change 
resulted in a small positive value of ACp (equation 1). The 
calculated AC^^'' and A Cf^^ values were +51 and +7.5, 
respectively. 



Binding parameters of 5-methyl uridine monophosphate 
(rlMP) 

To understand the biological significance of the presence of 2'- 
OH in TMP, we synthesized rTMP (Figure SI). Titration of the 
PAZ domain with rTMP produced a remarkable exothermic 
signal in comparison with dTMP (Figure S4). The binding affinity 
of rTMP was in the micromolar range with a A„ value of 
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Figure 2. ITC profiles of the binding of nucleotides with the PAZ domain bound with various nucleotides. 

doi:10.1371/journal.pone.0094538.g002 
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Table 3. The thermodynamic constants obtained by isothermal titration calorimetry for the association of nucleosides with 
Ago2PAZ domain. 



Substrate 


n 




Kc 


AH 


AG 


ms 


Estimated Kb 










(kcal/mol) 


(kcal/mol) 


(kcal/mol) 




UMP 


2±0.1 


6.6 ±0.6 


15 


-0.97±0.1 


-6.4 


5.3±0.3 




AMP-UMP 


2.7±0.2 


4.9 ±0.5 


20 


-0.24±0.2 


-6.5 


6.3±0.3 


K^"'' = U3 


dAMP-UMP 


2.7±0.2 


6.2 ±0.5 


62 


-0.68±0.5 


-6.6 


6±0.3 




dCMP-UMP 


2.1 ±0.1 


5.6±0.4 


17.8 


-0.75±0.6 


-6.6 


5.8±0.2 




Uridine-UMP 


1.2 ±0.02 


9.9 ±0.8 


10 


-1.9±0.05 


-6.9 


5.1 ±0.1 


K™'' = m 



The one set of sites binding model is used for fitting the data, (n = 3). 
doi:l 0.1 371 /journal.pone.C094538.t003 



17000 M indicating a 4-fold increase in binding affinity 
compared with dTMP. The interesting feature is the change of 
the enthalpic components of rTMP. While the binding of dTMP 
was enthalpically unfavorable (+520 cal/mol), rTMP binding 
showed favorable enthalpic conditions with a A// of — 64 cal/ mol. 

Gene silencing assay 

The gene silencing capability of siRNAs was estimated using a 
dual-luciferase assay. siRNAs that contained ribonucleotides 
showed more potent RNAi activity than siRNAs that carried 
deoxynucleotides at their 3 '-overhang (Fig. 3). Furthermore, the 
siRNA data bear an interesting correlation with data from ITC. 
The strongest RNAi was shown by rG or uridine, while low RNAi 
efficiency was produced by dC (Fig. 3). At 10 nM, all treatments 
showed a statistically significant RNAi effect compared with the 
control. Interestingly, at a 10-fold lower concentration (1 nM), 
only siRNAs bearing ribonucleotides showed statistically signifi- 
cant results compared with the control (Tables SI and S2). 

Discussion 

The PAZ domain is a universal binding pocket and adapts to a 
number of nucleotides with a wide range of binding aflinities 

The initial ITC experiments showed that some nucleotides 
(UMP, GMP, CMP, TMP, and dGMP) bound to the PAZ domain 
with considerable affmity. The lack of estimation of binding of 
some other ligands with the PAZ domain does not necessarily 
exclude the occurrence of a binding event. For instance, these 
substrates might bind with very low afiinity (beyond the accurate 
detection window of ITC); furthermore, the substrate might be 



negatively adsorbed onto the surface of the PAZ domain without a 
detectable exothermic interaction. To distinguish between these 
possibilities, the ITC experiments were repeated in the mutual 
displacement mode. At first, the binding constant and the 
thermodynamic parameters of the binding of the high affinity 
ligand (UMP) with the free protein were estimated. Then, we 
repeated the same experiment by titrating UMP in the presence of 
the PAZ domain bound with the low affinity ligand to check 
whether the affinity of UMP was changed in the presence of the 
other ligand. The difference between these two titrations can be 
used to detect the binding parameters of the low affinity ligand 
according the following equations [40]: K^pp = Ka/1+Kb[B] (3) 
and AH,pp = AHa-AHb(Ka/1+Kb[B]) (4), where Ka is the 
equilibrium constant for the tighdy binding ligand, Kg is the 
equilibrium constant for the weakly binding ligand, K^pp is the 
binding constant of the tightly binding ligand in the presence of 
the weakly binding ligand, A//^ is the enthalpy change accom- 
panying the weakly binding ligand and A//^ is the enthalpy change 
accompanying the tightly binding ligand and \E\ is the concen- 
tration of the free weakly binding ligand. 

According to the ITC-measured dissociation constants, the 
nucleotides binding with the PAZ domain can be classified into 
three grades. Grade 1 ligands: medium afTinity group includes 
UMP, GMP, CMP, and dGMP in which the affinity constant was 
> 1 0000 M ' . Grade 2 ligands: include low-medium affinity 
ligands such as TMP. Grade 3 ligands: are the low affinity ligands 
and include AMP, dAMP, and dCMP. 



Table 4. The thermodynamic constants obtained by isothermal titration calorimetry for the association of nucleosides with 
Ago2PAZ domain. 



substrate 


T 


Ka 




\H 


A5 


°C 


M^*10'' 


(HM) 


(cal/mol) 


(cal/mol) 


UMP 


8 


6.3 ±0.4 


15.8±0.5 


-1230±66 


21±1 




15 


1.6±0.2 


62.5 ±1.5 


-238±15 


18.4±1.1 




30 


1.3±0.1 


77±2.2 


-89±5 


12±1.2 


GMP 


8 


4.9±0.2 


20±1 


-251 ±20 


20.6 ±1 




15 


4±0.45 


25±1.2 


-173±10 


21.4±1.1 




30 


1.9±0.1 


52±4 


-84 ±0.5 


12.5 ±0.8 



The one set of sites binding model is used for fitting the data, (n = 3). 
doi:l 0.1 371/journal.pone.0094538.t004 
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Figure 3. Dual luciferase assay. The materials and conditions are described in the Methods section. *Significant at 0.05 level. **Significant at 0.01 
level. 

doi:1 0.1 371 /journal.pone.0094538.g003 



Mechanism of interaction of nucleotides with PAZ 
domain 

To elucidate the possible forces contributing to the recognition 
of nucleotides, we analyzed and correlated the binding site data 
with the thermal data. Alignment of the sequences of PAZ domain 
from different Argonaute proteins is given in Figure S5. The cavity 
enclosing the nucleotides (amino-acids enclosed in green boxes) is 
almost composed of hydrophobic residues with a minor contribu- 
tion from hydrophilic amino-acids. Furtheremore, the hydropho- 
bic residues constitute the direct contact with nucleotides (F68, 




Figure 4. Surface representation of the interaction of nucleo- 
tides with PAZ domain of Drosophila. Accessible surface 
representation of the Drosophila Ago2 PAZ domain showing the 
binding cavity of the Ago2 PAZ domain bound with CMP at its 3'- 
terminal (A). Ribbon representation of the PAZ domain bound with 
siRNAis shown in (B). Hydrogen bonds (shown in yellow) between the 
ribose of CMP and the side chain of GLN548 (C). The coordinates of 
interaction were derived from PDB ID 4F3T. 
doi:1 0.1 371/journal.pone.0094538.g004 



F86, HI 10, Till and Yl 12). Interestingly, the high conservation 
of nucleotide-binding amino-acids in PAZ domains indicates a 
common PAZ domain binding pattern among different Agos. 

The base or sugar of nucleotides, which is critical for 
binding with the PAZ domain? 

With respect to the size of the nucleotide, there was no 
discrimination based on their size. The small bases, such as 
cytidine, uridine, and thymidine in rCMP and UMP, were bound 
to the cavity of the PAZ domain with an affinity comparable to 
that produced by the larger purines, such as guanosine of GMP, 
Therefore, the size of the base attached to the nucleoside is not the 
critical factor for binding with the PAZ domain. On the other 
hand, we observe a different phenomenon when we compared the 
functional type of the sugar ring attached to the nucleic acid base. 
It was apparent that the PAZ domain has different binding 
preferences, such as being selective for ribonucleotides over 
deoxyribonucleotides. For instance, the affmities of dGMP and 
dCMP were far lower than those of their respective ribonucleo- 
tides. 

From the thermodynamic experiments, it was apparent that the 
absence of 2' -OH adversely affected the binding of nucleotides 
with the PAZ domain, probably due to the loss of important 
hydrogen bonding with the PAZ domain. This was manifested by 
the unfavorable enthalpic conditions accompanying dGMP and 
dTMP binding. The crystal structure of human EIF2C2/Ago2 
was resolved when either uridine or guanosine nucleotides were 
bound to the PAZ domain (PDB ID 4EI1 and 4F3T, respectively). 
Examination of the structure of the Drosophila Ago2 PAZ domain 
with ohgonucleotides reveals that the first nucleotide and, to a 
lesser extent, the second one, interact with the PAZ domain (Fig. 
4A), while the rest of the sequence faces the solvent (Fig. 4B). 
Within this configuration, the 2'-OH of the ribose ring interacts 
with glutamine (GLN548) through hydrogen bonding (Fig. 4C). 
Examination of the human EIF2C2/Ago2 PAZ domain structure 
bound with uridine nucleotide reveals that the 3 ' UMP is engulfed 
within a deep cavity of the PAZ domain (Fig. 5A). UMP adopts a 
U-shaped structure with the sugar ring at the bottom and is 
swallowed by the cavity. Furthermore, both the 2'- and 3' -OH 
groups of the sugar ring are bound to the PAZ domain through 
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Figure 5. Surface representation of the interaction of nucleotides with PAZ domain of Human Ago2. Accessible surface representation 
of the human EIF2C2/Ago2 PAZ domain (A) showing the binding cavity of the PAZ domain bound with UMP at its 3'-terminal. Stereoview showing 
the hydrogen bonding of the ribose sugar with the side chains of HIS1 1 0 and TYR1 1 2 (B). The binding cavity of the EIF2C2/Ago2 PAZ domain bound 
with GMP at its 3'-terminal (C). Hydrogen bonds between the ribose of GMP and the side chains of HIS1 10 and TYR1 12 are shown in green (D). The 
coordinates of interaction were derived from PDB ID 4EI1 and 4F3T. The figure was generated using Molsoft. 
doi:1 0.1 371 /journal.pone.0094538.g005 



bidentate hydrogen bonds with the side chains of HIS 1 1 0 and 
TYRI 12 (Fig. 5B). A similar configTaration was also found for the 
binding of guanosine nucleotides (Fig. 5C) in which bidentate 
hydrogen bonds were evident between the hydroxyl groups of the 
ribose and the side chains of HISllO and TYRI 12. Such 
interactions of the 2' -OH of ribose provide the basis for RNA 
selection and lower affinity of the PAZ domain for deoxyribonu- 
cleotides. 

In order to check our hypothesis that the nucleotide-bound 
sugar ring is the critical factor for EIF2C2/Ago2 PAZ domain 
binding regardless of the type of attached nucleic acid base, we 
synthesized rTMP and compared its binding parameters with 
dTMP. rTMP is an analogue of dTMP that contains a ribose 



sugar instead of the deoxyribose of dTMP. Thermodynamic data 
indicated drastic changes in binding affinity after the addition of 
2'-OH of the synthesized rTMP. Furthermore, there were 
favorable enthalpic conditions that probably arose from the ability 
of rTMP to make hydrogen bonds with the side chains of HIS 1 1 0 
and TYRI 12. The large difference in the enthalpic components of 
dTMP and rTMP suggests that other favorable enthalpic forces 
are more important than those formed between 2' -OH and the 
side chains of HISllO and TYRI 12. The presence of 2' -OH in 
ribonucleotides did not only facilitate binding through hydrogen 
bonds with the PAZ domain but also allowed a favorable 
conformation of the nucleotide for other binding sites, which are 
eventually lost with deoxyribonucleotides. Recently, the impact of 
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2' -OH modifications of siRNA on yeast Ago I activity was 
investigated [41]. Interestingly, tlie substitution of 2'-OH with 
2'-H at several positions, including the nucleotides binding to the 
PAZ domain, resulted in a significant impairment of yeast Agol 
activity. Together, these further clarify the role of 2'-OH not only 
for recognition by the PAZ domain but also for the overall activity 
of Ago proteins. 

Bias of the EIF2C2/Ago2 PAZ domain against adenine 
nucleotides 

The binding data revealed that all forms of adenine nucleotides 
showed the lowest binding efficiency among the tested nucleotides. 
The bias of cellular proteins against certain nucleotides was 

reported previously [42]. Specially referring to EIF2C2/Ago2, 
several reports indicated an unusual association between EIF2C2/ 
Ago2 and adenine nucleotides. In one study, the trimming 
behavior of a set of miR-124 mutant constructs with the 4 
different nucleotides at their 3 '-terminal was investigated [43]. 
Cells overexpressing EIF2Cl/Agol, EIF2C3/Ago3, and 
EIF2C4/Ago4 did not show inherent cleavage affinity to any of 
the terminal nucleotides. In contrast, EIF2C2/Ago2-overexpress- 
ing cells showed high cleavage activity of the mutant miR-124 that 
contained adenosine nucleotides at its 3 '-terminal. Furthermore, 
the cleaved product was one nucleotide shorter than the miR-124 
mutants containing G or C or U at their 3 '-terminal. Additionally, 
trimming of the 3'-terminal was in the order A<C<G<U. 
Consistent with the model provided by Juvvuna et al. [43], in 
which the highest exposure to trimming exonucleases in Ago2- 
loaded miRNA was associated with A nucleotides and the lowest 
rate was with U nucleotides, we assume that the very low affinity of 
A nucleotides allowed a weak interaction with the PAZ domain 
and they became vulnerable to exonucleases. In this regards, based 
on our thermodynamic data, we suggest a model for the cleavage 
of nucleotides at the 3'-terminal (Fig. 6). The presence of 
adenosine results in loose binding with the PAZ domain, making 
it vulnerable to cleavage by cellular nucleases. In contrast, the 
strong binding of U/G/C nucleotides does not allow the cleavage 
process to occur. In another study, the correlation between the 
adenylation and uridylation of RNAs with Ago-specific binding 
was investigated [44]. The authors reported that adenylation 
appeared to reduce the effectiveness of the miRNA targeting of 
mRNA transcripts. Furtiiermore, miRNA bound with EIF2C2/ 
Ago2 had a low number of adenine nucleotides at its 3'-terminal 
relative to miRNA extracted from the whole cell. Therefore, on 
the basis of our data, we suggest that there may be an usual 
terminal adenylation rate occurring at the 3 '-end of RNAs in the 
cell; however, the fraction bound with EIF2C2/Ago2 is the 
smallest among Ago proteins due to its inherent low affinity for 
adenine nucleotides. 

Combining these previous two examples together, either 
trimming or extension of the 3 '-terminal shared a common 
characteristic, that is the presence of an adenine nucleotide at this 
site infers a common fate. The presence of adenine nucleotides at 
the 3'-end of an miRNA improves its trimming or decreases its 
binding with the EIF2C2/Ago2 PAZ domain. 

Reviewing the previous observations as well as our ITC data, 
which indicate the lowest binding affinity for adenine nucleotides, 
there could be an inherent bias for the EIF2C2/Ago2 PAZ 
domain against adenine nucleotides. 

On the basis of the previously discussed data, we suggest a 
model outlining the possible interaction of nucleotides at the 3'- 
terminal of siRNA (Fig. 6). There are two possibilities for 
nucleotides at the 3 '-terminal of siRNA; either strong binders 
such as G, C, or U nucleotides (Fig. 6 Al) or weak binders such as 



an A nucleotide (Fig. 6 Bl). In the presence of strong binders, the 
PAZ domain and EIF2C2/Ago2 adopts a favorable conformation 
for proper RNAi (Al and A2). In contrast, the A nucleotide (B2) 
binds loosely or not at all to the PAZ domain, thereby becoming 
vulnerable to cleavage by nucleases and thus exposing the second 
nucleotide in the siRNA sequence (B3, B4). The process of 
cleavage can be repeated if the following nucleotide in the 
sequence is also an A nucleotide (CI and C2); otherwise, stronger 
binders can interact with the PAZ domain and induce RNAi (B5 
and B6). To this end, we would expect different and interesting 
interactions of siRNAs with the PAZ domain according to the 
nature of the nucleotides present at the 3 '-terminal of the siRNA. 

The EIF2C2/Ago2 PAZ domain is an induced-fit 
nonflexible binding cavity 

The energetic parameters of substrate-receptor binding allows 
us to estimate whether the binding mechanism operates by an 
induced fit process or is accompanied by conformational changes. 

The characteristic behavior of the induced fit process comes from 
the presence of high entropic barriers. The changes in AS along 
the measured experimental temperatures were very small, that is 
we were unable to measure Ts (the temperature at which the 
entropy of the system reaches zero). Therefore, empirical 
equations to derive the forces contributing entropy cannot be 
apphed in this case. Since entropy is the measure of the disorder of 
the system, a small change in entropy infers that no or small 
changes in conformation are observed in the experimental 
temperature range. In this context, major conformational changes 
are usually accompanied by large and positive entropic changes. 

The change in heat capacity is an indication of structural 
changes as it is associated with changes in tlu; surface area upon 
ligand binding. The change in heat capacity during the binding of 
GMP and UMP with the PAZ domain was very small and carried 
a positive sign, indicating that there were conformational changes 
accompanying substrate binding by the PAZ domain. Taken 
together, the small values of AS and ACp indicated that induced fit 
is the general mechanism for the binding of nucleotides with the 
PAZ domain. In order to vaUdate this assumption, we calculated 
the average changes in the surface components of the PAZ domain 
using the available structures (Table S3). The observed changes in 
surface area between the Apo and cytidine nucleotide-bounded 
Drosophila Ago2 PAZ domain were very small (+350 A). In 
addition, applying the polar and non-polar components to 
equation 2 yields a measured ACp of +6.4. Such a small calculated 
ACp value indicates the lack of significant conformational 
changes. A similar finding was observed after comparing the 
binding of the human EIF2C2/Ago2 PAZ domain with either U 
or G nucleotides, in which there was little or no change in the 
presence of various substrates. Such results correlate with the 
hypothesis that the RNA-binding cavity of the Agol PAZ domain 
is a preformed pocket 32. Thus, this feature seems to be conserved 
among Ago proteins. 

On the basis of these observations, we assume that during 
RNAi, the PAZ domain moves as a single unit or whole molecule 
toward the 3 '-end of the RNA with no or little internal 
conformational changes. 

Correlation between the measured binding affinity and 
the observed RNAi 

The obtained RNAi data indicated the general preference of 

RNAi for ribonucleotides. A surprising issue in these data was the 
in VIVO activity of siRNAs containing A nucleotides at their 3'- 
terminal, in spite of poor binding with the PAZ domain. In this 
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Figure 6. Model outlining the possible interactions of nucleotides at the 3'-terminal of siRNA. 

doi:10.1371/journal.pone.0094538.g006 
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context, siRNAs containing U or G were among the strongest 
inducers of RNAi at 1 and 1 0 nM, respectively. Within the known 
Ago proteins, the most catalytically active member is Ago2, while 
Agol, Ago3, and Ago4 are inactive. Ago2 removes the passenger 
strand promptly; however, a major problem in RNAi experiments 
is that the other Ago proteins remain bound with the double- 
stranded siRNA for several days, thus, they slowly accumulate the 
guide strand [45]. Furthermore, little is known about the 
functional properties of other Ago proteins and this might 
complicate the in vivo results. Additionally, blunt-ended siRNA 
duplexes are also capable of inducing a gene silencing effect 
[46,47]. When the role of the PAZ domain in RNAi was analyzed, 
it was reportedly essential for the specific and productive 
incorporation of siRNA into the RNAi pathway [21-23]. 
Importantly, PAZ-disrupted Ago mutants are unable to unwind 
or eject the passenger strand, thus highlighting the importance of 
the PAZ domain in the formation of a functional RISC complex 
[25]. Since ITC involves the direct measure of the association 
between two binding partners without any other interfering factor, 
we rely on ITC data to classifying the substrates according to their 
affmities. The effect of the strength of binding of siRNA with PAZ 
domain can be regarded as stage-specific matter. Stable complexes 
of 3' terminal of siRNA with PAZ domain are required during 
binding of dsRNA, during unwinding and release of the passenger 
strand. This is also important for guiding the RISC complex to 
bind with the target mRNA. In contrast, the release of 3' terminal 
of siRNA from PAZ domain is important for base pairing with 
target RNA and target cleavage. 

In conclusion, the EIF2C2/Ago2 PAZ domain showed 
selectivity for ribonucleotides over deoxyribonucleotides. The 
highest binding affinity was shown by UMP, CMP, and GMP. 
There was little change in entropic conditions and heat capacity, 
indicating a more or less static preformed cavity for siRNA 
binding. The PAZ domain showed an inherent bias against 
adenine nucleotides, which showed low binding affinity. The 
strongest RNAi efficacy was shown for siRNAs containing uridine 
or guanosine nucleotides at their 3'-terminal. Finally, conserva- 
tion of the 2'-OH or introducing a similar interaction seems to be 
essential for maintaining their high binding affinity with the PAZ 
domain. 



References 

1. Montgomery MK (2004) RNA interference: historical ovcr\'icw and significance. 
Mctliocls in molecular biology 265: 3—21. 

2. Wang QC, Nie QH, Feng ZH (2003) RNA interference: antiviral weapon and 
beyond. World journal of gastroenterology: VVJG 9: 1657—1661. 

3. Jinek M, Doudna JA (2009) A llrree-dimensional view of the molecular 
machinery of RNA interference. Nature 457: 405-412. 

4. Song JJ, Smith SK, Harmon GJ, Joshua-Tor L (2004) Crystal structure of 
Argonaute and its implications for RISC slicer activity. Science 305: 1434—1437. 

5. Riedmann LT, Schwentner R (2010) miRNA, siRNA, piRNA and argonautes: 
news in small matters. RNA biology 7: 133—139. 

6. tianey SA (2007) Expanding the repertoire of RNA interference screens for 
developing new anticancer drug targets. Expert opinion on therapeutic targets 
11: 1429-1441. 

7. Mathupala SP, Mittal S, Guthikonda M, Sloan AE (2007) MicroRNA and brain 
tumors: a cause and a cure? DNA and cell biology 26: 301—310. 

8. Zhang M, Wang J, Yao R, Wang L (2013) Small Interfering RNA (siRNA)- 
Mediated Silencing of the M2 Subunit of Ribonucleotide Reductase: A Novel 
Therapeutic Strategy in Ovarian Cancer. International journal of gynecological 
cancer: official journal of the International Gynecological Cancer Society. 

9. Bagasra O (2005) RNAi as an antiviral therapy. Expert opinion on biological 
tirerapy 5: 1463-1474. 

10. WooddeU CI, Rozcma DB, Hossbach M.John M, Hamilton HL, ct al. (2013) 
Hepatocyte-targeted RNAi Therapeutics for the Treatment of Chronic Hepatitis 



Supporting Information 

Figure SI Structure and NMR of 5-methyl uridine monophos- 
phate. 
(DOC) 

Figure S2 Thermodynamic signatures of UMP, GMP, CMP, 

TMP and dGMP binding with Ago2PAZ domain. 

(DOC) 

Figure S3 the binding aflfmity of UMP, GMP, CMP, TMP and 

dGMP binding with Ago2PAZ domain. 

(DOC) 

Figure S4 Titration of Ago2PAZ domain with the rTMP. The 
top panel shows the raw calorimetric data referring to the amount 
of heat produced following each injection. The bottom panel 

shows the integrated amount of heat generated per injection as a 

function of the molar ratio of rTMP to protein. 

(DOC) 

Figure S5 Amino acid alignment from different Agos. The 
residues forming the cavity of nucleotides binding are included in 
green boxes. The residues of direct contact with nucleotides and 
enclosed by red boxes. 
(DOC) 

Table SI The output of analyzed RNAi (1 nM) data by 
STATA. A single-tailed one-way analysis of variance with 
Bonferroni's multiple comparison test was conducted. 
(DOC) 

Table S2 The output of analyzed RNAi (10 nM) data by 
STATA. A single-tailed one-way analysis of variance with 
Bonferroni's multiple comparison test was conducted. 

(DOC) 

Table S3 The calculated accessible surface area parameters of 
Ago2PAZ domain. VADAR program is used in the calculations 
with a solvent probe radius of 1.4 A. 
(DOC) 

Author Contributions 

Conceived and designed the experiments: MK Y. Kitade. Performed the 
experiments: MK AA RN TS AK YN Y. Kilamura Y. Kitade. Analyzed 
the data: MK AA AK Y. Kitade. Contributed reagents/materials/anafysis 
tools: MK Y. Kitade. Wrote the paper: MK Y. Kitade. 



B Virus Infection. Molecular therapy: the journal of the American Society' of 
(jcne Therapy. 

11. Bian Z, Xiao A, Cao M, Liu M, Liu S, et al. (2012) Anti-HBV efficacy of 
combined siRNAs targeting viral gene and heat shock cognate 70. Virology 
journal 9: 275. 

12. Yang YJ, Zhao PS, Zhang T, Wang HL, Liang HR, et al. (2012) SmaU 
interfering RNAs targeting the rabies virus nucleoprotein gene. Virus research 
169: 169-174. 

13. Wu W, Wang T, DongJJ, Liao ZL, Wen FQ,(2012) Silencing of c-kit with small 
interference RNA attenuates inllammation in a murine model of allergic asthma. 
International journal of molecular medicine 30: 63—68. 

14. Popescu FD, Popescu F (2007) A review of antisense therapeutic interventions 
for molecular biological targets in asthma. Biologies: targets & therapy 1: 271- 
283. 

15. Courties G, PresumeyJ, Duroux-Richard I, Jorgensen C, Apparailly F (2009) 
RNA interference-based gene therapy for successful treatment of rheumatoid 
arthritis. Expert opinion on biological therapy 9: 535—538. 

16. Kubowicz P, Zelaszczyk D, Pekala E (2013) RNAi in Clinical Studies. Current 
medicinal chemistry. 

17. Ramachandran PV, Ignacimuthu S (2013) RNA Interference-A Silent but an 
Efficient Therapeutic Tool. Applied biochemistry and biotechnology. 

18. Sasaki T, Shiohama A, Minoshima 8, Shimizu N (2003) Identification of eight 
members of the Argonaute family in the human genome small star, filled. 
Genomics 82: 323-330. 



PLOS ONE I www.plosone.org 



12 



May 2014 I Volume 9 | Issue 5 | e94538 



EIF2C2/Ago2 PAZ Domain Interactions 



19. Song JJ, Joshua-Tor L (2006) Argonautc and RNA — getting into the groove. 
Current opinion in structural biology 16: 5—11. 

20. lingcl A, Simon B, Izaurralde E, Sattler M (2004) Nucleic acid 3'-end 
recognition by the Argonaute2 PAZ domain. Nature structural & molecular 
biology 11: 576-577. 

21. Ma JB, Ye K, Patel DJ (2004) Structural basis for overhang-specific small 
interfering RNA recognition by the PAZ domain. Nature 429: 3 1 8-1^22. 

22. SongJJ, LiuJ, Tolia Ml, Schneidernian J, Smilh SK, et al. (2003) The crystal 
structure of the Argonaiite2 PAZ domain reveals an RNA liinding motif in 
RNAi eflcctor complexes. Nature structural liiology 10: 1026 1032. 

23. Yan KS, Yan S, Parooq A, Han A, Zcng L, et al. (2003) Structure and conserved 
RNA binding of the PAZ domain. Nature 426: 468-474. 

24. Lingel A, Simon B, Izaurralde E, Sattler M (2003) Structure and nucleic-acid 
binding of the Drosophila Argonautc 2 PAZ domain. Nature 426: 465—469. 

25. C;u S, Jin Huang Y, Zhang P, Kay MA (2012) SlicingTndependent RISC 
Activation Requires the Argonautc PAZ Domain. (.Current biolog}': CB 22: 
1536-1542. 

26. Taniho K, Nakashima R, Kandecl M, Kitamura Y, Kitadc Y (2012) Synthesis 
and biological properties of chemically modified siRNAs bearing 1-deoxy-D- 
ribofuranose in their 3'-overhang region. Bioorganic & medicinal chemistry 
letters 22: 2518-2521. 

27. Ueno Y, Kawada K, Naito T, Shibata A, Yoshikawa K, et al. (2008) Syndiesis 
and silencing properties of siRNAs possessing lipophilic groups at their 3'- 
termini. Bioorganic & medicinal chemistry 16: 7698—7704. 

28. Somoza A, Terrazas M, Eritja R (2010) Modified siRNAs for the study of the 
PAZ domain. Chemical communications 46: 4270-4272. 

29. Liu J, Carmell MA, Rivas FV, Marsdcn CO, Thomson JM, et al. (2004) 
Argonaute2 is the catalytic engine of mammEilian RNAi. Science 305: 1437— 
1441. 

30. Maiti M, Nauwelaerts K, Lescrinier E, Herdewijn P (2011) Structural and 
binding study of modified siRNAs with the Argonautc 2 PAZ domain by NMR 
spectroscopy. Chemistry 17: 1519-1528. 

31. Zeng L, Zhang Yan K, Zhou MM (2011) Structural insights into piRNA 
recognition by the human PlWMike 1 PAZ domain. Proteins 79: 2004-2009. 

32. Tian Y, Simanshu DK, MaJB, Patcl DJ (201 1) Structural basis fbr piRNA 2'-0- 
methylated 3 '-end recognition by Piwi PAZ (Piwi/Argonaute/ZwiUe) domains. 
Proceedings of the National Academy of Sciences of the United States of 
America 108: 903-910. 

33. Simon B, Kirkpatrick JP, Eckhardt S, Renter M, Rocha EA, et al. (2011) 
Recognition of 2'-0-methylated 3'-end of piRNA by the PAZ domain of a Piwi 
protein. Structure 19: 172-180. 



34. Kandeel M, Kitadc Y (2013) Computational Analysis of siRNA Recognition by 
the Ago2 PAZ Domain and Identification of the Determinants of RNA-Induced 
Gene Silencing. PloS one 8: e57140. 

35. Kandeel M, Kitade Y (2013) In sihco molecular docking analysis of the human 
Argonautc 2 PAZ domain reveals insights into RNA interference. Journal of 
computer-aided molecular design 27: 605-614. 

36. Kandeel M, Kitade Y (2011) Binding dynamics and energetic insight into the 
molecular forces driving nucleotide binding by guanylate kinase. Journal of 
molecular recognition: JMR 21: 322 332. 

37. Kandeel M, Kitadc Y ;2010) Sulistrate speeilicil\ and nucleotides liinding 
properties of NM23H2/nucleoside diphosphate kinase homolog from Plasmo- 
dium falciparum. Journal of bioencrgeties and biomembranes 42: 361-369. 

38. Kandeel M, Kitade Y (201 1) The substrate binding preferences of Plasmodium 
thymidylate kinase. Biological & pharmaceutical bulletin 34: 173-176. 

39. Kandeel M, Miyamoto T, Kitade Y (2009) Bioinformatics, cnzymologic 
properties, and comprehensive tracking of Plasmodium falciparum nucleoside 
diphosphate kinase. BiologicEd & pharmaceutical bulletin 32: 1321-1327. 

40. Sigurskjold BW (2000) Exact smalysis of competition Hgand binding by 
displacement isothermal titration calorimetry. Analytical biochemistry 277: 
260-266. 

41. Nakanishi K, Weinberg DE, Bartcl DP, Patel DJ (2012) Structtire of yeast 
Argonautc with guide RNA. Nature 486: 368-374. 

42. Ho T, Wang L, Huang L, Li Z, Pallett DW, ct al. (2010) Nucleotide bias of DCL 
and AGO in plant anti-virus gene silencing. Protein & cell 1: 847—858. 

43. Juwuna PK, Khandelia P, Lee LM, Makeyev EV (2012) Argonautc identity 
defines the length of mature mammalian microRNAs. Nucleic acids research 40: 
6808-6820. 

44. Burroughs AM, Ando Y, de Hoon MJ, Tomaru Y, Nishibu T, et al. (2010) A 
comprehensive survey of 3' animal miRNA modification events and a possible 
role for 3' adcnylation in modulating miRNA targeting effectiveness. Genome 

research 20: 1398 1410. 

45. Petri S, Ducck A, Lchmann Piitz N, Rude! S, el al. (201 1) Increased siRNA 
duplex stability correlates with reduced ofl-target and elevated on-target eticcts. 
RNA 17: 737-749. 

46. Ucno Y, Inouc T, Yoshida M, Yoshikawa K, Shibata A, et al. (2008) Synthesis of 
nuclease-resistant siRNAs possessing benzene-phosphate backbones in their 3'- 
overhang regions. Bioorganic & medicinal chemistry letters 18: 5194-5196. 

47. Ueno Y, Watanabe Y, Shibata A, Yoshikawa K, Takano T, et al. (2009) 
Synthesis of nuclease-resistant siRNAs possessing universal overhangs. Bioor- 
ganic & medicinal chemistry 17: 1974-1981. 



PLOS ONE I www.plosone.org 



13 



May 2014 I Volume 9 | Issue 5 | e94538 



